Ground-state properties of bosons in three- and two-dimensional traps 



Augusto Gonzalez 1 ' 2 !^ and Aurora Perez 2 ! 

1 Departamento de Fisica, Universidad de Antioquia, AA 1226, Medellin, Colombia 
2 Institute de Cibernetica, Matematica y Fisica Calle E 309, Vedado, Habana 4, Cuba 

(February 1, 2008) 



OO 
G\ 

CO 



> 

CO 
CN 

o 

CN 
O 
OO 

o\ 



c 

o 
o 



X 



We study trapped systems of bosons at zero temperature 
in three and two dimensions. Conditions are fulfilled for the 
application of Gross-Pitaevskii theory with a positive scatter- 
ing length. Series expansions for ground-state properties are 
obtained in both the noninteracting and the strong-coupling 
(Thomas-Fermi) limits. From these expansions, analytic es- 
timates are presented in the form of two-point Pade approx- 
imants. We explicitly show the approximants for the total 
energy per particle and the chemical potential. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj 



Since the discovery of Bose-Einstein condensation in 
alkali-vapour atomsa, trapped bosonic systems have at- 
tracted a lot of attention. In the experiments, the num- 
ber of confined atoms, N, ranges between 10 4 and 10 6 , 
whereas the ratio between the scattering length and the 
harmonic oscillator length is a/ a^o ~ 10~ 3 , thus the con- 
ditiona-fpr the application of the Gross-Pitaevskii (GP) 
theorytm are fulfilled. We assume that the interaction 
between pairs of atoms is repulsive, i.e. a > 0, as for 
87 Rb. The temperature is taken to be zero, that is we 
are dealing with ground-state properties. 

The GP equation for the condensate function, iji, in 
three dimensions is written as 



2m 



A + V ext (r)+g\^\ 2 (f)-fi}iP = 0, 



(1) 



where m is the mass of the atoms, fi — the chemical poten- 
tial, and g = 4-7T h 2 a/m. The condensate function satisfies 
the constraint 



d 3 r IV'I 2 = N. 



(2) 



In (1-2), we have neglected any effect coming from par- 
ticles out of the condensate. V ex t(r) is the external po- 
tential responsible for the confinement of the atoms. For 
simplicity, we will study an isotropic trap, that is V ext = 
^muj 2 r 2 , and is a symmetric function ip(r) — ip(r). 

A scaling of variables, r — > a ho r, — > (N/a^) 1 ^ 2 ^, 
fx — » fiujfi, reduces eqs. (1,2) to the dimensionless form 



-A+-r 2 +g\ip\ 2 -Li\ip = 



d 3 r |0| 2 = 1, 



(3) 



(4) 



and makes explicit that any scaled magnitude will de- 
pend only on the variable g = Ng/ (huja^) — AirNa/aho, 



where a^o = h 1 ' 2 /( 
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This scaling is preserved to 
the extent the GP equation remains valid. For example, 
the u> can not be increased up to a value at which a^ 
becomes comparable to a. 

In terms of the condensate function, the chemical po- 
tential is written as 
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d d rW--A + -<h/> + ffM 



(5) 



It differs from the total energy per particle in half the 
Hartree energy 



E= / d 



1 



1 



9 1 



0*( A+-r^ + ^l 



(6) 



We will consider formally that the variable g ranges be- 
tween zero (noninteracting bosons) and infinite (Thomas- 
Fermi theory). In the <? ^ limit, we may apply pertur- 
bation theory, i.e. to look for ij) and /i in the form 

= 00 + "013 + • ■ ■ > (?) 
fi = fx + nig + (8) 

The leading contributions are given byi@ 
7r -3/4 e -r /2 ^ Q _ 3^2. Next corrections are easily ob- 
tained also, resulting in 
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(n|0 2 |O>, 



n>0 
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Mi - (0|^o 2 |0) 



(27T) 3 / 2 ' 



(9) 
(10) 



where the \n) are three-dimensional harmonic oscillator 
states. 

The series for E is obtained from (6): Eq = /UOj whereas 
Ej. = mi/2. 

On the other hand, in the formal g — > oo limit, the ki- 
netic energy may be neglected and the so-called Thomas- 
Fermi (TF) theory for bosons appliestra. The conden- 
sate function is given by -0oo = sj (R 2 — r 2 )/(2g), where 
R = (155/(47r)) 1 / 5 , and the chemical potential reads 
Moo = R /2. The leading contribution to the energy is 
Eoo = 5R 2 /U. 

The boundary layer near the condensate surface is re- 
sponsible for the corrections to the TF theoryO, leading 
to contributions of the order of 1/R 2 , 



1 



E = h R2 + ^ Hlm2R) 



2^1n(1.413 R) 



(11) 
(12) 



Taking together the weak-coupling and strong- 
coupling series, for anyone of the magnitudes E and /i 
we can write 



f(R)\ 



b + b 5 R 5 + 0(R 10 ), (13) 
R 2 {a + ^HAR) + O(l/R 5 )}. (14) 



The coefficients bo, 65, ao and 04 are listed in Table 1. 
Notice that the variable R coincides with the conden- 
sate radius of the TF theory only in the g — > 00 limit 
(differently from the notation usecLin Ref. [5]). 

In most experimental conditions^], interactions are nei- 
ther so weak to be considered in perturbation theory, nor 
so strong for the TF theory to be valid. In the present 
paper, two-point Pade approximants {P s j{R)} are to be 
constructed as analytic estimates to the magnitudes f(R) 
over the entire range of variation of R. s + 1 coefficients 
from the expansion (20) and t + 1 from (21) are used 
to determine the coefficients in the approximant. No- 
tice that in (13-14), we have at our disposal in total 15 
coefficients (many of which are equal to zero). 

By construction, the approximants are asymptotically 
exact in both the R — > and R — » 00 limits. We will show 
that the error of the higher approximants is lower than a 
few percents at any R. A recent application of the Pade 
technique to obtain the ground-state energy of electrons 
in a parabolic quantum dot led to similar resultsDll 

As it is usual in Pade-approximant techniqueso, we 
shall show convergence of a sequence {P s j}, in which 
both s and t increase simultaneously. In the present prob- 
lem, the sequence {Pk+3,k} exhibits good convergence 
properties. The explicit expressions for the first nontriv- 
ial terms of this sequence are the following 



P 5)2 (P) = b 



hR 5 



1 + qiR + q 3 R 3 ' 
qi = b b 5 /al, <?3 = b 5 /a , 

5 1 + qiR 



P 6 , 3 (R) = bo + b 5 R b 



l + qi R+... + q 4 R 4 
q-3 = b 5 /a , qi = 3q 3 /(2a ), 
94 = hqi/a , q 2 = 3b 5 qi/(2al). 



(15) 



(16) 



We show in Fig. 1 the relative difference between P$ t2 
and Pq i3 for the total energy per particle (the K = 2 
curve). The maximum relative error of Pq j3 may be esti- 
mated from this curve to be less than 7% over the entire 
range of variation of R. Analogous results are obtained 
for the chemical potential. 

To build up the next approximant, we need the ~ 1/ R 2 
term in the high-P series. However, this term contains 
a logarithmic function of R. We may circumvent this 
problemE3 by constructing a (7,4) interpolant for the 
magnitude 
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76 5 i? 5 



OiR 1 



(17) 



R — >oc 



4a n P 2 



R 2 



0{1/R 3 ). (18) 



The series (17-18) are similar to (13-14) up to a re- 
definition of coefficients: bo = 2&o, ^5 = 76s, ao = 4ao, 
a4 = 04. The approximant takes the form 



1 + Q 1 R + Q 2 R 2 



(19) 



P 7A (R) = b + b 5 R 5 

1 + QiR + . . . + Q 5 R J 
where Q 5 = b 5 Q 2 /a , Qi = a Q 2 /feo, 
Q3 = (ao/Q 2 + bo)Qs/ao, 

Qi = QiQ 2 / Q5, 

and Q 2 is obtained from 

a Qi - boQ 3 + = 0. 

The magnitude f(R) is obtained by multiplying by R 
and integrating the resulting expression . We will call it 
again the approximant, 



Pta(R) = b 



7&5 
R 2 



dx x 6 



1 + Qtx + Q 2 x 2 



l + Qix- 



Q 5 x 5 



(20) 



The integration over x could be explicitly performecO, 
but a direct numerical integration of (20) is trivial. The 
coefficients Qi, . . . , Q5 are listed in Table 1. 

The relative difference between P74 and P@ 3 is also 
shown in Fig. 1 (the K = 3 curve). The maximum error 
of P7 4 may thus be estimated to be < 1.2%. A similar 
result is obtained for /1. 

Next, we turn to the two-dimensional situation. We 
may think of an anisotropic three-dimensional trap in 
lu 2 z 2 )/2, where r 



which V e . 
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and u! z /u! Xt y >> 1, so that the motion of bosons in 
the z-direction is described by a gaussian of very small 
width. Such—highly anisotropic traps have been already 
constructed^. 

In eq. (1), we write ty(r,z) = x( z )' t P( r )i where x( z ) = 



l/4„-r 



/ ( 2h ) Multiplying the equation by 



X and integrating over z, the resulting equation takes 
again the form (1), but the parameters entering it are 
9 = 93Djdz x i = 93D y / muj z /(2TT), fi = fi 3D - ^w 2 /2, 
and all the integrations will run over two-dimensional 
space. The reference level for any other magnitude with 
dimensions of energy will be huj z /2 also. 

A scaling of variables, r — > dh r, ip — ► (iV/a^) 1 / 2 ^, 
/i — > Tuj/j,, in which ui = oj x>y , reduces the GP equation 
to the dimensionless form (3), where g = Ng/(huja 2 ). 
Notice that, differently from the 3D case, g does not de- 
pend on ijj. Thus, the dependence on u> of any magnitude 
can be obtained on purely dimensional grounds. For ex- 
ample, any energy is exactly proportional to Tilj. This is, 
of course, valid to the extent the GP equation is valid. 
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In the g limit, the solution is looked for as (7- 
leading to 



V^o=^ 1/2 e- r2/2 , Mo = l, 
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n>0 
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Mi - WolO) 



2tt 



(21) 

(22) 
(23) 



Now, the \n) are two-dimensional harmonic oscillator 
states. The coefficients of the series for E are, again, 
E Q = /i , Ei = Hi/2. 

On the other hand, in the g — > 00 limit, TF theory 
leads to Voo = y/{R 2 -r 2 )/(2<?), where i? = (Ag/n) 1 ^ 4 , 
/.too = R 2 /2, and = i? 2 /3. Corrections are again 
given by boundary layer theory. We performed in 2D 
calculations similar to those of Ref. [5] . The results are 



1 9 
E = -R 2 

O 

1 9 
= -i? 2 



3i? 2 



ln(1.604 i?) 
ln(1.435 R) 



(24) 
(25) 



2" ' 3i? 2 

For anyone of the magnitudes E and /i, we have then 
f(R)\ R ^ = b + b 4 R 4 + O(R s ), (26) 
l^oo = R 2 {«o + ^HAR) + 0(1/R 5 )} . (27) 

The coefficients are listed in Table 2. 

Pade approximants are to be constructed from (26-27). 
Notice that, in the approximation we are working, f(R) 
contains only even powers of R or 1/R. Thus, the max- 
imal power in the Pade should be even. In the sequence 
{Pk+5,k}, this leaves the approximants 



P 5 ,o{R) = b o + 



b 4 R 4 



Pr,2(R) = bo + hR 



1 + (b 4 /a )R2 ' 

1 + q 2 R 2 



(28) 



1 + q 2 R 2 + q±R 4 
q 2 = Mo/(ao - 64), 94 = b 4 q 2 /a Q . (29) 

The relative difference between J^o and P?2 for the 
energy per particle in two dimensions is shown in Fig. 2 
(the K = curve). The maximum relative error of P? t 2 
may be estimated to be < 9% (< 10% for /i). 

The next approximant in this sequence, the P9.4, makes 
use not only of the a 4 coefficient, but of the &g as well, 



P%a{R) = b a + 



664 
1$ 



dx 



I + Q2X 2 



1 + Q2X 2 + Q 4 x 4 + Q 6 x 6 ' 



(30) 



where P 4 = {ao/b4)Qe, and Q2, Q4, Qe are obtained from 
the equations 



64Q4 + b s = a Q 6 , 
(I0Q4 — boQa = 64Q2, 
aoQ 2 — boQ 4 ~b 4 — a 4 Qe, 



(31) 
(32) 
(33) 



The modified coefficients are 60 = 26o, b 4 = 664, bs = 
lO&s, ao = 4ao, a 4 = a 4 . The numerical values for Q 2 , 
Q 4 , Qe, and P 4 are listed in Table 2. 

The coefficient 6s is obtained from second order per- 
turbation theory For the chemical potential the result 



H2 = 3(^0^1 



3 V WoW 

9. n 



(34) 



n>0 



The coefficient E 2 is obtained from (6), which may be 
rewritten in the g — » limit as 



E 



Mi - 2/i 2 ^ 2 



(35) 



The relative difference between P9 4 and P-j t2 is drawn 
also in Fig. 2. It is labelled as the K = 2 curve. From 
this curve, we may estimate the maximum relative error 
of the Pg,4 approximant to be less than 1.8% at any R. 

A direct comparison with the numerical calculations of 
Ref. [12] for the energy per particle in two dimensions in 
the thermodynamic (N — > 00) limit is presented in Fig. 
3. The difference between both results is always below 
the predicted 1.8%. 

In conclusion, we obtained analytic Pade estima- 
tions to ground-state properties of bosons in the Gross- 
Pitaevskii theory. The approximants work with a small 
relative error (< 2%) at any boson density. As examples, 
we explicitly found the approximants for the energy per 
particle and the chemical potential in both three and two 
dimensions. 

Pade approximants for other magnitudes can be con- 
structed in the same way. Different trap geometries could 
be considered also. Finite systems, for which the GP for- 
malism is no longer valid, could be studied also because 
we can easily apply perturbation theory near g = 0, 
and use a Hartree approximation in the g — > 00 limit. 
Finite temperatures could be equally well considered, 
in particulac, because of the scaling of thermodynamic 
magnitudesEJ. All these analytic estimations could be of 
great value for the experimental groups. 



ACKNOWLEDGMENTS 

The authors acknowledge support from the Colombian 
Institute for Science and Technology (COLCIENCIAS) 
and from the Research Grant Programme of the Third 
World Academy of Sciences (TWAS). Both authors are 
associate members of the International Center for Theo- 
retical Physics (ICTP, Trieste). 



3 



* Electronic mail: agonzale@democritus.udea.edu.co 
' Electronic mail: aperez@cidet.icmf.inf.cu 
1 M. H. Anderson, J. R. Ensher, M. R. Mathews, C. E. Wie- 

man, and E. A. Cornell, Science 269, 198 (1995); K. B. 

Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. 

S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 

75, 3969 (1995); C. C. Bradley, C. A. Sackett, J. J. Tollett, 

and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 
2 E. P. Gross, Nuovo Cimento 20, 454 (1961); L. P. 

Pitaevskii, Sov. Phys. - JETP 13, 451 (1961). 

3 S. Giorgini, L. P. Pitaevskii and S. Stringari, J. Low Temp. 
Phys. 109, 309 (1997). 

4 F. Dalfovo and S. Stringari, P hys. Rev. A 53, 247 7 (1996). 

5 A. L. Fetter and D. L. Feder, |cond-mat/9704173 . 

6 M. J. Holland, D. S. Jin, M. L. Chiofalo, and J. Cooper, 
Phys. Rev. Lett. 78, 3801 (1997). 

7 A. Gonzalez, J. Phys.: Condens. Matter 9, 4643 (1997). 

8 A. Gonzalez, B. Partoens, and F. M. Peeters, Phys. Rev. 
B 56, 15740 (1997). 

9 G. A. Baker, Jr. and P. Graves-Morris, Pade approximants 
(Addison-Wesley, Reading, Mass., 1981). 

10 A. Ishihara and E. W. Montroll, Proc. Nat. Acad. Sci. 68, 
3111 (1971). 

11 M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. M. 
Kurn, D. S. Durfee, C. G. Townsend, and W. Ketterle, 
Phys. Rev. Lett. 77, 988 (1996); C. J. Myatt, E. A. Burt, 
R. W. Christ, E. A. Cornell, and C. E. Wieman, Phys. Rev. 
Lett. 78, 586 (1997). 

12 T. Haugset and H. Haugerud, 

13 S. Giorgini, L. P. Pitaevskii and S. Stringari, Phys. Rev. 
Lett 78, 3987 (1997). 



TABLE II. Coefficients in 2D. 



:ond-mat/9707107 









bo 


1 


1 


bi 


1/16 


1/8 


b s 


-0.001124 


-0.003371 


Q,0 


1/3 


1/2 


CL4 


4/3 


2/3 


Q2 


0.646107 


0.751621 


Q 4 


0.292455 


0.424002 


Qe 


0.073825 


0.142145 


Pi 


0.262488 


0.379052 



FIG. 1. Relative differences between consecutive approxi- 
mants for the energy in three dimensions. 



FIG. 2. The same as in Fig. 1, but in two dimensions. 



FIG. 3. Comparison between the Pg t 4 approximant and the 
numerical calculations of Ref. 12 for the energy per particle 
in two dimensions. 



TABLE I. Coefficients in 3D. 





E 


M 


bo 


3/2 


3/2 


b 5 


l/(15\/27r) 


2/(15v / 2tt) 


ao 


5/14 


1/2 


0,4 


5/2 


3/2 


Qi 


0.302359 


0.30288 


Q 2 


0.082748 


0.084582 


Q 3 


0.152967 


0.209793 


Qi 


0.039404 


0.056388 


Qs 


0.010784 


0.015747 
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